Equilibration of a one-dimensional Wigner crystal 
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Equilibration of a one-dimensional system of interacting electrons requires processes that change 
the numbers of left- and right-moving particles. At low temperatures such processes are strongly 
suppressed, resulting in slow relaxation towards equilibrium. We study this phenomenon in the case 
of spinless electrons with strong long-range repulsion, when the electrons form a one-dimensional 
Wigner crystal. We find the relaxation rate by accounting for the Umklapp scattering of phonons 
in the crystal. For the integrable model of particles with inverse-square repulsion, the relaxation 
rate vanishes. 

PACS numbers: 71.10.Pm 



The low energy properties of systems of interacting 
fermions in one dimension are commonly described in the 
framework of the Luttinger liquid theory [l|, Q • This the- 
ory successfully predicted a number of interesting proper- 
ties of one-dimensional electron systems, such the power- 
law renormalization of the tunneling density of states and 
impurity potential 0, 0]. On the other hand, much re- 
cent interest in interacting one-dimensional Fermi sys- 
tems was focused on the phenomena not captured by 
the Luttinger liquid theory [^-Q. One example involves 
equilibration of a moving one-dimensional electron liq- 
uid, which was recently shown to affect the conductance 
of quantum wires Q and drag between two wires ■ 

In the case of weakly interacting one-dimensional elec- 
trons the physical mechanism of equilibration was dis- 
cussed in Ref. ll|. At low temperature T excitations 



of the system are particle-hole pairs near the two Fermi 
points, with the typical energy T <^ Ep and momentum 
T/vp ^ Pf- (Here Ep, vp, and pp are the Fermi energy, 
velocity and momentum of the system.) In this regime, 
backscattering of a right-moving electron near the Fermi 
level requires transfer of momentum 2pp to a large num- 
ber of particle-hole pairs. The most efficient such process 
consists of a sequence of scattering events, in which a hole 
passes from the left to the right Fermi point through the 
bottom of the band. Such processes are suppressed as 
^-Ep/T g^jj^fj^ consequently, equilibration of the chemical 
potentials of the right- and left-moving electrons is a very 
slow process. 

Equilibration of one-dimensional fermions beyond the 
weak interaction regime is a more challenging problem. 
In this case the description in the language of particles 
and holes is no longer applicable, and according to the 
Luttinger liquid theory the elementary excitations of the 
system are bosons. On the other hand, the Luttinger 
liquid theory does not adequately describe particles near 
the bottom of the band, and is therefore incapable of 
describing the equilibration processes. In this paper we 
show that this difficulty can be overcome in the limit of 



strong long-range interactions. 

More specifically, we consider a system of identical 
spinless particles of mass m described by the Hamilto- 
nian of a general form 



I 1,1' 



(1) 



Here xi and pi are the coordinate and momentum of the 
l-ih particle, and V{x) is the interaction potential. In 
the limit of very strong repulsion, the particles form a 
periodic chain with interparticle distance a — \/n deter- 
mined by their density n. In the case of Coulomb repul- 
sion, V{x) = /\x\, such an arrangement is commonly 
referred to as the Wigner crystal. 

At strong but finite repulsion, the particles can deviate 
from their respective lattice sites, xi = la + ui, but the 
relative change of interparticle distance remains small, 
\ui — ui>\ <C \l — l'\a. To leading order in the deviations 
ui the Hamiltonian ([TJ takes the form 



(2) 



{Ui ~ Uvf 



(2) 



l.V 



where we used the following notation for the r-th deriva- 
tive of the interaction potential 



V, 



(r) _ 



(FV{x) 



dx'' 



(3) 



x—la 



Elementary excitations of the harmonic chain ^ are 
phonons characterized by quasimomcntum q (i.e., ui oc 
^iqiy Their frequencies are easily found by solving the 
classical equations of motion. 



^l/(^)[l-cos(gO]. 



(4) 



Provided that the interaction potential V{x) falls off 
faster than l/\x\ at large distances, the low-energy exci- 
tations of the system are bosons with acoustic spectrum 
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FIG. 1: Umklapp scattering of two right-moving phonons 
(filled circles) into two left-moving phonons (open circles). 
The energies of the initial and final states are equal (see the 
dashed circle). On the other hand, the sum of the quasimo- 
menta defined to range in the first Brillouin zone — tt < q < tt, 
decreases by 27r. 



u}„ = s\q\, in agreement with the Luttinger liquid theory 



Importantly, in the limit of strong interactions the har- 
monic approximation ([2]) provides the full spectrum of 
elementary excitations (j?]), not limited by the restriction 
\q\ <^ 1 imposed in the Luttinger liquid theory. An- 
other advantage of the Wigner crystal model is that the 
weak interaction of phonons is naturally described the 
anharmonic terms in the expansion of the Hamiltonian 
([T]) in the powers of the displacements ui. Scattering of 
phonons resulting from these interactions leads to relax- 
ation of their distribution function to equilibrium. 

At temperatures T much lower than the Debye energy 
fkjjTr , the quasimomenta of the thermally excited phonons 
are small, q 1. For such phonons the Umklapp scatter- 
ing is impossible, and phonon-phonon collisions conserve 
the total quasimomentum Q of the system. As a result, 
the equilibrium distribution of the phonons 



^1 ^h{ui,-uq)/T _ I 



(5) 



is characterized by two parameters: the temperature T 
and the velocity u of t he phonon gas with respect to the 
lattice (see, e.g., Ref. 



13|) 



It is important to note that while being a good approx- 
imation, conservation of quasimomentum is not exact. 
Indeed, even at low temperatures there is a finite occu- 
pation iV^ ~ g-fiw^/T phonon states near the edge 
g = TT of the Brillouin zone, leading to a small probabil- 
ity of Umklapp processes. As a result, the total quasimo- 
mentum of the phonons relaxes as Q = —Q/t with the 
time constant t c>c e^^^l'^ , For distribution ([5]) one has 
Q (X u, and thus the velocity u acquires a time depen- 
dence, ii = —u/t. Calculation of the relaxation time t is 
our main goal. 

The microscopic mechanism of Umklapp scattering is 
illustrated in Fig. [TJ We assume that the phonon spec- 
trum is concave, which is the case for Coulomb repul- 



sion. Then the dominant process involves scattering of a 
rare phonon with quasimomentum qi near the boundary 
of the Brillouin zone by an acoustic phonon with energy 
Hujq^ ^ T. As a result of such a collision the quasimomen- 
tum of the high-energy phonon changes by Sqi ~ T/hs. 
If the new quasimomentum q'l = qi + 6qi is outside the 
Brillouin zone (— tt, tt), the scattering involves Umklapp, 
and the total quasimomentum changes by ±27r. 

In the course of such scattering events the number of 
the rare phonons near the edge of the Brillouin zone is 
conserved, while their momentum changes by a small 
amount <C tt. Thus these phonons essentially diffuse 
in the momentum space, and the evolution of their distri- 
bution function P{q,t) can be described by the Fokker- 
Planck equation 



dtP = -d„J, 



J = A{q)P{q,t) ~ -d,[B{q)P{q,t)]. (6) 



Here J has the meaning of the probability current in 
momentum space, and the functions A{q) and B{q) can 
be expressed in terms of the rate IV^^g+s^ of the phonon 
transition from state q to q + Sq as 



A{q) = ^ SqWg,g+Sq, B{q) - ^(<5<z)^M^,,,+5, 

6q Sq 



(7) 



Below we will be using the Fokker-Planck equation ([5]) to 
determine the behavior of the phonon distribution func- 
tion P in a small vicinity of q = tt. This enables us 
to approximate B(q) hy B = B(tt). Furthermore, we 
will be interested in the case of weak deviation from 
equilibrium, when the velocity of the phonon system is 
small, u ^ s. Thus we can approximate A{q) by its 
value at equilibrium, when the Boltzmann distribution 
P{q) — e^''"''/^ has to solve the Fokker-Planck equation 
by nullifying the probability current J. Then, from ([6]) 
we find A{q) = —Bfvjj'^l2T , where = dujq/dq. 

The phonon distribution function P{q) is 27r periodic. 
In the first Brillouin zone and away from its edges q — ±7r 
it is given by Eq. ([S]). A periodic continuation of ([5]) 
results in a formal discontinuity at q = tt. Specifically, 
the phonon distribution function away from q — t: takes 
the form 



P{q) = e-'i'^',/T^±TThu/T^ 



T 



< T{q - tt) < TT. 



(8) 

We now solve the Fokker-Planck equation ([6|) with the 
boundary conditions (|S]) to determine the steady-state 
phonon distribution P{q) at \q - n] < yr/^^- The 
solution is rather straightforward and analogous to that 
for the distribution function of holes near the bottom of 
the conduction band of weakly interacting electrons 11 1 . 



It interpolates smoothly between the Boltzmann func- 
tions (|H]) and corresponds to a finite but exponentially 
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small probability current 

The non-vanishing value of J means that in unit time 
{N/2'k)J phonons increase their momentum and leave 
the first Brillouin zone through point q — t:. (Here 
N is the total number of electrons in the Wigner crys- 
tal.) Each such event is an Umklapp process resulting 
in the decrease of the quasimomentum Q of the phonon 
system by 27r. We therefore conclude that at non-zero 
phonon velocity u the phonon scattering events result in 
Q — —NJ. Comparing this result with the total quasi- 
momentum of the phonons Q — ttNuT'^ /3s'^ , easily com- 
puted using the distribution ([5|), we find the relaxation 
rate = —Q/Q in the form 

As expected, the relaxation rate shows activated temper- 
ature dependence, with activation temperature given by 
the Debye energy of the phonons huj-j^. However, to fully 
determine the temperature dependence of the prefactor, 
one has to calculate the diffusion coefficient B of phonons 
in momentum space. 



It is convenient to treat the scattering of phonons using 
the second quantization of the Hamiltonian ([1]) whereby 
the displacements and momenta of the particles are ex- 
pressed in terms of the phonon destruction and creation 
operators bq and as 

The quadratic Hamiltonian ([5]) then takes the standard 
form 

Ho = E ^^9(^9^9 + 1/2). (13) 

The coupling of phonons is described by the anharmonic 
corrections to Hq, which are easily obtained by expanding 
the full Hamiltonian ([1]) to higher orders in the displace- 
ments u/. To find the leading contribution to phonon 
scattering, we will need to account only for the cubic 
and quartic anharmonisms, H « Hq + U^^^ + U^^\ with 
the respective perturbations taking forms 



t/(3) 



3/2 



/a (91, 92) 



q^^q^ V^gi'^<Z2 '^91+92 

/4(gl,g2,g3) 
48m2iV ^^^^^ ^UJq,UJq^UJq^Uq,+q^+q^ 



-92 + ^gi+92)' 



(14) 



E 



(V + blqjibq, + blqjibq, + blqj (b. q, .q, .q, + b^+q^+qj. (15) 



Here the functions /s and are defined as 



/algi, ?2) = ^i ^H^Mili + q2)l] - sin((ji/) - sin(q20}: 

(16) 



1=1 



and 



/4(gi, 92,93) = ^ V^^''^{1 - cos((ji/) - COS((720 
1=1 

- cos((j3/) - cos[(gi +q2+ q3)l] 

+ cos[(gi + q2)l] + cos[(gi + gs)/] 

+ C0s[{q2 +q3)l]}. (17) 

The rate of two-phonon scattering processes shown in 
Fig. [TJ in which a phonon qi moves to the state q'l, is 



given by the golden rule expression 



9l,9i 



27r 



E K^02^l\.,'fNqANq'^+l) 
92,92 



92->9i,92l ^'92V^'92 
^^qi+q2,q[+q'2^i^qi +^92 -^91 -^q^' (1^) 



The scattering of two phonons qi and q2 into q'l and q'2 , 
can be accomplished in the first order in the quartic an- 
harmonism U^'^\ Alternatively, the same scattering pro- 
cess can be realized in second order in the cubic anhar- 
monism U^^\ Simple power counting shows that in both 
cases the resulting amplitude is proportional to h?, i.e., 
one has to account for both contributions. The actual 
calculation is straightforward and results in the scatter- 
ing matrix element in the form 



t. 



A 



9l,92^9i:9i 



m^N [iOq^UJq^UJq'^UJq'^Yl'^ ' 



(19) 
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where 



A 



qi+q2 



92- 91 



92 ^qi , 



■'9i> 



To find the dffiusion constant B{qi) of phonons in mo- 
mentum space, one notes that the momentum (72 is hm- 
ited by the occupation number Nq^ in Eq. (jlSp , resulting 
in \q2\ ~ |(?2l ^ T/hs <^ tt. Thus one can expand A in 
powers of Sq — q[ — qi = 92 — 92- The expansion starts 
with a quadratic term, A cx (Sq)'^. The proportionality 
constant depends on the specific model of the interaction 
potential V{x). In the most interesting case of Coulomb 
potential V{x) — e^/|x| our treatment is complicated 
by the fact that the phonon speed s diverges logarith- 
mically. In practice, however, the Coulomb potential is 
always screened at large distances by remote gates. In 
this case, A can be found analytically. For simplicity, we 
also limit ourselves to the most important case of gi = tt 
and find 



(21) 



where (^(x) is the Riemann's zeta function. Combining 
this result with Eqs. (HH), (HI]), and ©, we find 



B = xT^ 



21^3^(3) 



20 



(22) 



One can now substitute this result into Eq. (ITOl) to 
obtain the full expression for the relaxation rate of 
the phonon system in a Wigner crystal. In the 
specific case of Coulomb interaction, the Debye fre- 
quency ujj, — (7C(3)e^/ma3)^/^, parameter uj'^ — 
— 2 ln2(e^/7C(3)ma'^)^/^, and the speed of phonons is 
s — {2e'^\n{d/a)/ma'^y/^, where d is the distance to the 
gate. 

Although the expression (|22l) is derived for Coulomb 
interaction, the temperature dependence B oc is valid 
for any long-range repulsive potential. An interesting 
exception is the case of strong inverse-square repulsion, 
V{x) = with 7 ^ f? jm. This form of repulsive 

potential corresponds to the Calogero-Sutherland model, 
which is exactly solvable due to the presence of an infinite 
number of integrals of motion. As a result, one expects 
that scattering of excitations preserves their momenta, 
and no diffusion in momentum space should be possible. 
Indeed, we have been able to show that the coefficient 
X in Eq. ([22t vanishes for inverse-square repulsion. We 
have also verified that the expression ((20|) for A, and thus 
the phonon scattering amplitude (|19p. vanish in this case. 
Scattering of a massive particle off an acoustic phonon. 
Fig. m was recently discussed in the context of quantum 
decay of dark solitons in one dimensional Bose systems 



14| . In analogy with our observation, their decay rate 
vanishes in the integrable (Lieb-Liniger) case. 

The expression (fTO|) for the relaxation rate of phonon 
system in a one-dimensional Wigner crystal with the dif- 
fusion constant ((22|) is the main result of this paper. To 
illustrate its significance we now briefiy discuss the ef- 
fect of the phonon relaxation upon the conductance of 
a strongly-interacting quantum wire. A more rigorous 
treatment will be presented elsewhere. 

The equilibration of phonons in a one-dimensional 
Wigner crystal is analogous to electron equilibration in 
the limit of weak interactions. In the latter case, the ex- 
citations are holes/particles created by transferring elec- 
trons to/from a Fermi point. The momentum of each ex- 
citation is thus measured from the nearest Fermi point. 
Equilibration processes include backscattering of holes 
near the bottom of the band. Since the momenta of right- 
and left-moving holes are measured from different Fermi 
points, each backscattering event changes the momen- 
tum of the excitations by 2pi? — lix ja^ in analogy with 
the quasimomentum change AQ — 27r when a phonon in 
a Wigner crystal crosses the edge of the Brillouin zone. 
We thus conclude that Q can be interpreted as the rate of 
backscattering of right-moving electrons, = —Q/2tt. 

This relation provides for a way to observe the phonon 
equilibration by measuring conductance of a strongly in- 
teracting wire connected to non-interacting leads. Neg- 
ative value of means that some of the right-moving 
electrons entering the wire from the left return to the 
same lead, thereby reducing the conductance [gI. [qI. [ii| . 
Adapting the calculation of Ref. 11] to the case strong 
interactions, one finds 



G = 



1 - L 



B /h^\L 



1/2 



(23) 



Although the correction to the conductance is exponen- 
tially small, it grows with the length of the wire L. The 



result ((23)) applies to relatively short wires, 
^huj^/T^ c.f. Ref. [Hi. 



cx 



To summarise, we have studied the relaxation of the 
distribution function of phonons in a one-dimensional 
spinless Wigner crystal. Full equilibration of phonons 
requires Umklapp processes, resulting in an exponen- 
tially small relaxation rate. The preexponential factor 
in the resulting expression PU)) scales as T^/^. In the in- 
tegrable case of inverse-square interactions the relaxation 
rate vanishes. 
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